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Abstract. We present near-infrared narrow- and broad-band imaging of the field of GRB 021004, performed with 
ISAAC on the UTl of the ESO Very Large Telescope. The narrow-band filters were chosen to match prominent 
emission lines at the redshift of the absorption-line systems found against the early-time afterglow of GRB 021004: 
[O III] at 2=1.38 and Ha at 2=1.60, respectively. For the 2=1.38 system we find an emission-line source at an impact 
parameter of 16", which is somewhat larger than the typical impact parameters of a sample of Mgll absorbers at 
redshifts around unity. Assuming that this tentative redshift-identification is correct, the star formation rate of 
the galaxy is 13±2 M© yr~^. Our study reaches star-formation rate limits {5a) of 5.7 Mq yr^^ at z=1.38, and 
7.7 Mq yr~^ at 2=1.60. These limits correspond to a depth of roughly 0.13 L*. Any galaxy counterpart of the 
absorbers nearer to the line of sight either has to be fainter than this limit or not be an emission-line source. 

Key words, gamma rays: bursts - galaxies: distances and redshifts - quasars: absorption lines 



1. Introduction 

Identification of the objects responsible for the high 
redshift intervening absorption systems in QSO spec- 
tra, such as Damped Lya (DLA) and Mgii absorbers, 
has been the expressed goal of a large number of 



observing prograinmes over the past 2 decades (e.g. 
Smith et al. 1989; "Dcharveng et al.' "1990'; jLe Brun etall 



1997t iGuillemin fc Bergeron, .1997t .Warren et alJ l200ll 



Kulkarni et alJ I2OOII: rColbert fc MalkanI 1200^ Even 
though many such programmes have been successful in 
detecting objects at (or near) the absorption redshift, a 
doubt does in almost all cases persist as to whether the 
identified object is indeed the absorber or if there could be 
another object hidden inside the glare of the QSO itself 
at much smaller impact parameter. It was only realized 
very late that this problem is seriously exacerbated by 
the strong bias of the absorption cross-section selection 
which picks out mostly object s with the smallest impact 
param eters llFvnbo et alJll999l:ISchavell200lUMoller et all 
l2002bt) . The discovery of bright, high-redshift optical tran- 
sients (OTs) of gamma-ray bursts (GRBs) has now made 
it possible to completely eliminate this last doubt. GRB 
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OTs, like QSOs, allow identification of high-redshift ab- 
sorbers, but the OTs, unlike QSOs, fade away completely 
within a few months leaving the spectroscopically investi- 
gated sightline open for extremely deep imaging studies. 

The current status concerning the galaxy counter- 
parts is the following. For the high-z (z>2) DLA sample, 
the debate has b een between advocates for large, disk- 
formi ng galaxies ijWolfe et alJ Il986t IProchaska fc Wol^ 
200lh and faint and small, gas-rich dwarfs l)Tvsonlll98a 
Haehnelt e t al.'l998V n three cases DLA absorbing galax- 
ies have been identified as faint blue Lya emitting galaxies 
at impact parameters ranging from 0.99 to 2.51 arcsecs (8- 
19 kpc), supporting the view that DLA galaxies at z>2 
are faint blue dwarf galaxies ()Moller et al. 2002b,b In the 
range from z=0 to z=1.2, Guillemin & Bergeron (1997) 
reported that galaxy counterparts of Mg 11 absorbers with 
equivalent widths larger than 0.6 A predominantly are spi- 
rals with luminosities similar to present-day L* galaxies. 
The impact parameters of the identified galaxies were 
found to correlate with the luminosities of the galaxies, 
such that bright galaxies are found at larger impact pa- 
rameters. In a flux limited study, objects below the given 
detection limit will fail to be identified, and because of 
the correlation with impact parameter this will tend to 
bias all identifications both towards brighter objects and 
towards larger impact parameters. An indication that this 
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Fig. 1. The transmission curves for the four filters used in this 
project. From left to right it is Js, NB 1.19, H, NB 1.71. 



may be the case for the Mgll studies is the clear tendency 
for higher redshift absorbers to be identified as brighter 
galaxies with larger impact parameters (around 110 kpc 
at z=l.l) than those at lower redshifts (around 20 kpc at 
z=0.1). The actual absorbers may therefore still be hiding 
at small impact parameters below the detection limit. 

GRB 021004 is the third brightest afterglow de- 
tected so far, with R =15.4 only 4 minutes after the 
burst iIFox et al.ll2003l l. Spectroscopic observations (e.g. 
IChornock &: Filippen ko 2002) found a host galaxy redshift 
of z=2.33. In addition two other systems were identified 
with absorption lines of Mgi, Mgli (both systems with 
equivalent widht s larger than 0.6 A) and Fell at z=1.38 
and z=1.60 fe.g. iMaller et al.ll2002a|) . In this pilot study 
we present observations in narrow- and broad-band filters 
to search for emission lines from the galaxy counterparts 
of these absorbers. The emission-line imaging technique 
that we employ allows identification of the foreground ab- 
sorber, even when it would be located exactly along the 
line of sight to the GRB and its host galaxy. Throughout 
the paper we assume Hq = 70kms~^ Mpc~^, fim = 0.3 
and ilx = 0.7. 



2. Observations 

At z=1.38, [Oiii] 15007 is redshifted to 1.19 fim, and 
at z=1.60. Ha emission is shifted to 1.71 /im. These 
wavelengths correspond to the central wavelengths of two 
ISAAC narrow-band filters. Imaging in these two filters 
and in the Js and H band allows identification of the inter- 
vening systems if they are star-forming galaxies. In Fig.^ 
we show the filter transmission curves for the four filters. 

The observations were carried out in service mode with 
UTl/ISAAC of ESO's Very Large Telescope on the nights 
of December 17/18 and 19/20, 2002, for a total effec- 
tive exposure time of 51 and 105 minutes in the narrow- 
band filters NB 1.19 and NB 1.71, respectively. To enable 
continuum subtraction, we also observed the field in the 
broad-band filters Js and H, each for 11 minutes. Table 
m shows the log of the observations. The ISAAC short- 
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wavelength detector is a Rockwell 1024x1024 chip, with 
a gain of 4.5 electrons/ADU, a read noise of 11 electrons 
and a pixel-size of 0'.'148. 

From the raw images possible ghosts of especially 
brigh t objects were removed with the ECLIPSE l)Devillardl 
l200lh routine ghost. After dark subtraction, the flat-field 
for each image was constructed using the ECLIPSE rou- 
tine twflat. After flat-fielding and removal of bad pixels us- 
ing the fixpix task within IRAF^, the sky was determined 
for each image from the adjacent images on either side. In 
the case of the narrow-band filters, we used four images 
on either side, and for the broad-band filters we used as 
many as possible (10). We performed two sky passes. In 
the first pass all object images of a particular Observation 
Block (OB) are sky subtracted. These images are regis- 
tered with a shift only (i.e. no scaling nor rotation) to 
the coordinate system of the first image of the OB using 
about 5 bright objects in the field, and combined. From 
this image we make an object mask using the optiona l 
output image from SExtractor l|Bertin fc Arnout£lll996|) . 
which returns the pixels above the noise of the detected 
objects. With imreplace in IRAF, these objects were trans- 
formed into a bad pixel mask, making them larger by a few 
pixels. Then we perform another sky run, this time with 
all objects masked while average-combining the adjacent 
sky images. Again we register the images, and now com- 
bine them using the sigclip rejection in imcomhine within 
IRAF, which excludes possible cosmic rays and bad pixels 
that were missed. Finally, residuals left from bias varia- 
tions and possible horizontal jumps were removed by sub- 
tracting the median of each row and/or column of the 
image using ECLIPSE routines. This median was deter- 
mined from the final image with the objects masked. The 
final reduced images were trimmed to a size of 2' by 2', 
which is the part of the image that was exposed in each 
sub-exposure (the size of the original images is 2f5 by 2.'5, 
while we set the jitter box width to 30"). 

We use d observations of th e NICMOS standard 9123 
(S427-D) l)Persson et a,l.l l 



119981) to calibrate the narrow 
bands and observations of the UK IRT faint standards FS6 

to calibrate the Js and 



and FS32 (|Hawarclen 
H bands. For the broad bands, no correction for a colour 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Fig. 2. A 40x40 arcsec^ field centred on the position of GRB 021004 from the sum of the Js and H band images (left), the 
NB 1.19 image (middle), and the NB 1.71 images (right). The only emission-line source closer in than 30" is a candidate [O ill] 
emitter at an impact parameter of 16" seen near the western edge of the NB 1.19 image. 



term was attempted. We corrected for atmospheric extinc- 
tion assuming an extinction of 0.06 mag per unit airmass. 
We estimate that the uncertainties on the zero-points are 
0.04 and 0.10 for Js and H respectively. The zeropoints for 
the narrow-band filters are determined by integrating the 
spectrum of S427-D over the filter transmission curves to 
obtain the flux in erg cm^^ outside the atmosphere. 



3. Emission-line candidates 

Inspection of the images reveals no emission in any of 
the filters at the location of the early-time afterglow of 
GRB 021004. To select emission-line candi dates we fol- 
low a method similar to t he one outlin ed in 'Fvnbo et alJ 
l)l999(l . We let SExtractor ijBertin fc Arnouts detect 
objects using a very low threshold, i.e. requiring that 2 ad- 
jacent pixels are at least 1.5 sigma above the sky. For all 
detected sources we perform aperture photometry with 
the phot routine in IRAF, after which we reject all objects 
with a signal-to- noise ratio less than 5. The aperture is set 
to a radius of 6 pixels (0"9 on the sky) for all images. This 
roughly corresponds to 1.5 times the size of the full width 
at half maximum (FWHM) of the average image point- 
spread function (PSF). The local sky was determined in 
an annulus 10 pixels wide, at a distance of 10 times the 
image FWHM (about 40 pixels) from the object centre. 
The coordinates of the objects found on the narrow-band 
images were transformed to the broad-band images, and 
at these locations aperture photometry was performed in 
exactly the same manner, irrespective of an object being 
present or not. 

In both narrow-band images we detect a number of 
emission-line candidates (6 in NB 1.19 and 4 in NB 1.71 
above 5cr), but almost all at very large impact parameters 
(>30"). Only one candidate in the NB 1.19 image, i.e. 
[Oiii] at 2=1.38, has a small enough impact parameter to 
be a possible galaxy counterpart of the absorber consis- 
tent w ith the sample identified bv lGuillemin fc Bereeronl 
l)l997(l . It is located at 16" from the GRB line-of-sight and 



can be seen at the western edge of the section of the field 
shown in Fig. [S] In Fig. |21 we show the spectral energy 
distributions of the seven continuum objects named A-F 
in Fig. 121 and of the [O III] emission- line candidate named 
O Ilia. Although the emission-line nature of this object is 
clear, follow-up spectroscopy is needed to verify the red- 
shift. 

The NB 1.19 emission-line flux of O iiia is 8.1 x 10"" 
erg s~^ cm~^, while its continuum Str upper limit, esti- 
mated from the Js and H upper limits, is 2.1 x 10~^^ 
erg s~^ cm~^ in the NB1.19 fllter. Under the assumption 
that the emission line is indeed [Oiii], the candidate is 
at a luminosity distance of 3.0 x 10^^ cm and has an 
[Oiii] luminosity of Ljom] = 7 x lO''^ erg s~^. Using the 
er npirical flux ratio [O II I ] /Ha of 1/2 (based on Table 5 
of iGuillemin fc Bereeronl l)l997|) ). we obtain an estimated 
Ho; luminosity of 1.4 x 10"*^ erg s~^. This corr esponds to 
a sta r-formation rate (SFR) of 13±2 Mq yr~^ lKe n nicut3 
Il998|) . The error given here is the measurement error; 
an additional error of order 30% originates from the un- 
certainty in the Ha luminosity to SFR conversion factor 
l Kennicuttlll998>l . 

Our observations are sensitive to a limiting luminos- 
ity (5(J with an aperture radius of 6 pixels, and assuming 
the continuum contribution in the narrow-band filters is 
negligible compared to the emission line) of L[oin] = 3.6 

10*1 erg s-i 



10*1 erg s-i at z=1.38, and Lh^ = 9.7 



at 2=1.60. The Ha luminosity is about 0.13 L* relative to 
the Schechter parametrisation of the H a luminosi t y func - 
tion of galaxies at similar redshifts of lYan et alJ l|l999l) . 
These luminosity limits correspond to star-formation rates 



of SFR^ 



I Kennicut 



ti fil998^ 



5.7 Mo yr-i and SFR,=i.6o=7.7 Mq yr" 



4. Discussion 

The 16" impact parameter of the NB 1.19 candidate 
corresponds to 135 kpc at z=1.38. The sample of 
quasar absorption-line systems identified with galaxies by 
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iGuillemin fc Bergeron! l)l997|) were found to have separa- 
tions of 1.9-23" which, in the cosmology that we assume, 
correspond to 12-127 kpc. Our candidate is therefore con- 
sistent with, but in the upper range of, the impact param- 
eters of the Mgii systems previously studied. However, we 
find 9 emission-line objects down to a limiting SFR of 7.7 
Mq yr~^ in two fields each of angular size 2'x2', so there 
is a priori 50% chance that one of those objects would be 
found within an impact parameter of 16" in one of the 
fields. 

In co nclusion, this pilot study is c onsistent with the 
result of IGuillemin fc BergeronI l)l997|) that the typical 
Mgll absorber is an L* galaxy at an impact parameter 
of 100 kpc. In our case the candidate absorber is less lu- 
minous, but in the upper end of the impact parameter 
distribution. In their intermediate redshift (0.7-1.2) sam- 
ple, Guillemin & Bergeron identified 10 of 16 absorbers 
but failed to find counterparts to the remaining 6. We 
find one but fail to identify a candidate in the other field, 
while our star-formation limit (see above) compared to 
theirs (SFR^^o.T-i. 3=0. 7-43 Mg yr'^) is such that we 
would have expected to detect half of the objects they 
identified. 

The question therefore remains open whether these 
large galaxies are really the absorbers, or if they are merely 
the brightest neighbour to dwarf galaxy absorbers. For 
DLA absorbing galaxies at higher redshifts impact pa- 
rameters are found to be of order 10 kpc rather than 100 
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Fig. 3. Spectral energy distributions (propotional to Fa) of 
the seven objects A-F marked in Fig.|5|and for the candidate 
2=1.38 [Olll] emitter. The candidate, named Ollla and shown 
at the bottom, has strong excess flux in the NB 1.19 filter. 



kpc jMoller et al.ll2002b|) . and the SFRs of order 5-7 Mq 
yr~^. This is just below the 5a limit of our current data. 
A future larger study of GRB intervening Mg ii and DLA 
absorbers should aim at reaching at least a factor of three 
deeper (in SFR) in order to answer the question if we are 
still missing a population of faint absorbers at low impact 
parameters, and if Mgii absorbers are in fact related to 
dwarf galaxies of the DLA type. 
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